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LONG-TERM  GOALS 

The  long  term  objective  of  our  research  for  the  “High  Resolution  Air-Sea  Interaction”  (HIRES) 
Departmental  Research  Initiative  (DRI)  is  to  identify  the  couplings  between  large  wave  events, 
winds,  and  currents  in  the  surface  layer  of  the  marine  boundary  layers.  Turbulence  resolving 
large  eddy  simulations  (LESs)  and  direct  numerical  simulations  (DNSs)  of  the  marine  atmo¬ 
spheric  boundary  layer  (MABL)  in  the  presence  of  time  and  space  varying  wave  fields  will  be 
the  main  tools  used  to  elucidate  wind-wave-current  interactions.  A  suite  of  turbulence  simula¬ 
tions  over  realistic  seas  using  idealized  and  observed  pressure  gradients  will  be  carried  out  to 
compliment  the  field  observations  collected  in  moderate  to  high  winds.  The  database  of  simu¬ 
lations  will  be  used  to  generate  statistical  moments,  interrogated  for  coherent  structures,  and 
ultimately  used  to  compare  with  HIRES  observations. 


OBJECTIVES 

Our  near  term  goals  are:  1)  participate  in  the  planning  process  for  the  HIRES  field  campaign; 
and  2)  construct  an  LES  code  applicable  to  the  HIRES  high  wind  regime.  In  order  to  accom¬ 
plish  the  latter  goal  we  are  improving  the  parallelization  of  our  base  LES  code  and  developing 
an  algorithm  to  allow  simulations  of  turbulent  winds  over  nearly  arbitrary  3-D  wave  fields. 


APPROACH 

We  plan  on  investigating  interactions  between  the  MABL  and  the  connecting  air-sea  interface 
using  both  LES  and  DNS.  The  waves  will  be  externally  imposed:  (1)  based  on  well  established 
empirical  wave  spectra;  or  (2)  ultimately  provided  by  direct  observations  of  the  sea  surface  from 
field  campaigns.  The  main  technical  advance  is  the  development  of  a  computational  tool  that 
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allows  for  nearly  arbitrary  3-D  wave  fields,  i.e.,  the  sea  surface  elevation  h  =  h(x,y,t )  as  a  sur¬ 
face  boundary  condition.  The  computational  method  will  allow  time  and  space  varying  surface 
conditions  over  a  range  of  wave  scales  (P(10)m  or  larger. 


WORK  COMPLETED 

Meetings:  We  attended  PI  meetings  in  La  Jolla  CA  that  focused  on  refining  elements  of  the 
HIRES  held  campaign.  The  discussion  concentrated  on  the  atmospheric  measurements  to  be 
collected  from  aircraft  and  from  R/V  FLIP  with  particular  attention  paid  to  the  surface  layer 
pressure  measurements.  We  also  spent  several  days  at  the  Center  for  Interdisciplinary  Remotely- 
Piloted  Aircraft  Studies  (CIRPAS),  in  Monterey  CA,  while  the  HIRES  held  campaign  was  un¬ 
derway  in  June  2010. 

Algorithmic  Developments  and  Simulations:  During  the  past  year  we  completed  adapt¬ 
ing  our  parallel  LES  code  (Sullivan  &  Patton,  2008)  to  accommodate  a  general  3D  time  varying 
wavy  surface.  This  was  accomplished  in  the  following  stages:  a)  Develop  the  governing  equa¬ 
tions  for  an  atmospheric  PBL  in  curvilinear  time-dependent  surface-following  coordinates;  b) 
Build  and  test  a  module  capable  of  generating  a  3D  wavy  surface;  and  c)  Implement  and  test 
a)  and  b)  in  our  “hat  bottom”  parallel  LES  code  for  the  atmospheric  PBL.  Extensive  testing 
of  the  code  is  underway.  A  description  of  the  code  and  hrst  results  are  described  in  Sullivan  et 
al.  (2010a).  We  note  that  the  formulation  is  general  and  also  allows  turbulence  simulations  over 
and  around  3D  stationary  orography  as  described  in  Sullivan  et  al.  (2010b)  1 . 

Here  we  briehy  highlight  a  process  study  that  examines  the  influence  of  wave  age  on  the  dy¬ 
namics  of  the  marine  atmospheric  surface  layer.  In  these  simulations  the  geostrophic  wind  is 
varied  from  Ug  =  (5,  7.5, 10, 15,  20)  m  s_1  for  a  neutrally-stratified  marine  PBL  in  a  domain 
(. Xl,Yl,Zl )  =  (1200, 1200,  800)  m  using  (Nx,Ny,Nz)  =  (512,  512, 128)  gridpoints.  Thus  the 
horizontal  grid  spacing  Ax  =  Ay  =  2.34  m  and  the  hrst  vertical  level  is  1  m  above  the  water. 
The  synthetic  surface  waveheld  is  built  to  match  a  Pierson  and  Moskowitz  (1964)  wave  spec¬ 
trum  with  random  phases.  A  directional  spectrum  is  picked  to  emphasize  long  crested  waves. 

An  x  —  y  view  of  the  3D  surface  wave  held  is  given  in  hgure  1  and  an  x  —  z  slice  of  the  compu¬ 
tational  mesh  above  this  waveheld  is  given  in  hgure  2.  The  initial  temperature  sounding  6  =  300 
K  up  to  the  inversion  height  Zi  =  400  m,  beyond  this  height  6  increases  linearly  at  3  x  10~3  K 
m_1.  The  surface  heating  Q*  =  0,  the  surface  roughness  za  =  0.0002  m,  and  the  Coriolis  pa¬ 
rameter  /  =  10~4  s_1.  The  waveheld  is  built  based  on  the  assumption  of  a  surface  wind  speed 
of  15  m  s_1  and  the  phase  speed  of  the  peak  in  the  spectrum  Cp  ~  18  m  s_1.  Thus  the  suite 
of  simulations  allows  us  to  examine  a  wide  variation  of  wave  age  from  swell  dominated  to  near 
wind-wave  equilibrium.  Table  1  lists  bulk  properties  of  the  simulations,  viz.,  the  geostrophic 
wind,  wave  age,  and  friction  velocity  u*.  Dio  is  the  reference  wind  speed  at  a  height  of  10  m. 

The  simulations  are  run  for  more  than  50,000  timesteps  using  restart  volumes  with  fully  devel¬ 
oped  turbulence.  The  iteration  count  in  the  pressure  Poisson  solver  is  typically  set  to  30  and  the 
calculations  run  on  either  512  or  1024  computational  cores. 


1The  meeting  papers  Sullivan  et  al.,  (2010a, b)  are  available  at  http://www.mmm.ucar.edu/people/sullivan/. 


Table  1:  Simulation  properties 


Run 

Ug  (m  s  x) 

Cp/Uio 

u*  (nr  s  *) 

A 

5 

4.8 

0.124 

B 

7.5 

3.4 

0.187 

C 

10 

2.8 

0.228 

D 

15 

1.9 

0.338 

E 

20 

1.5 

0.452 

RESULTS 

Previous  field  observations,  turbulence  closure  modeling,  and  our  own  idealized  LES  (see  discus¬ 
sion  in  Sullivan  et  al. ,  2008)  all  show  that  fast  moving  swell  can  induce  marked  changes  in  the 
atmospheric  surface  layer  winds,  viz.,  an  upward  momentum  flux  from  the  ocean  to  the  atmo¬ 
sphere,  a  low-level  wind  maximum,  and  departures  from  law-of-the-wall  scaling.  The  prelimi¬ 
nary  LES  computations  performed  here  over  a  more  realistic  sea  surface  are  in  good  qualitative 
agreement  with  the  previous  studies  but  suggest  the  impact  of  swell  on  the  surface  layer  winds 
is  sensitive  to  the  content  of  the  wave  spectrum. 

One  of  the  surprising  results  from  the  present  simulations  is  the  significant  impact  of  swell  on 
the  coherence  and  magnitude  of  the  near-surface  pressure  fluctuations.  This  is  illustrated  in  fig¬ 
ure  3  where  we  compare  p'  /  p  for  two  levels  of  wind  forcing  Ug  =  (5,  20)  m  s-1,  i.e.,  a  low-wind 
situation  with  swell  and  a  high  wind  case  approaching  wind- wave  equilibrium.  The  difference 
in  the  pressure  signals  is  striking  and  even  more  remarkable  in  animations  of  the  pressure  held. 
In  the  low-wind  swell  case  there  is  a  very  strong  correlation  between  p' / p  <  0  and  wave  crests 
and  similarly  between  p' / p  >  0  and  wave  troughs  that  extends  over  the  depth  of  the  surface 
layer.  Inspection  of  the  how  visualization  and  animations  reveals  that  the  strong  correlation 
persists  across  the  range  of  resolved  waves,  i.e.,  both  large  and  small  scale  waves  appear  to  in¬ 
duce  a  similar  pressure  pattern.  The  coherence  of  the  wave  induced  pressure  held  can  extend  to 
20  m  or  more  depending  on  the  amplitude  of  the  underlying  wave.  Also,  the  pressure  signatures 
propagate  at  the  speed  of  the  waveheld,  additional  evidence  that  the  signals  are  generated  by 
surface  waves  and  not  atmospheric  processes.  These  are  clear  signatures  of  “wave  pumping”  by 
the  surface  waveheld  on  the  atmosphere.  The  amplitude  of  the  wave  spectrum  (and  hence  the 
level  of  wave  forcing)  is  held  constant  in  our  simulations  but  the  magnitude  of  the  turbulence, 
as  measured  by  w*,  increases  substantially  with  increasing  wind  speed.  The  structure  of  the  near 
surface  pressure  held  is  a  result  of  these  two  competing  effects.  At  low  winds  coherent  pressure 
signals  are  generated  by  the  wave  motions  when  the  turbulence  is  weak  but  this  coherence  is  de¬ 
stroyed  by  strong  turbulence  at  higher  winds. 

Figure  4  shows  that  the  impact  of  wave  age  also  appears  in  the  vertical  velocity  helds.  In  the 
low- wind  swell  regime  we  observe  large-amplitude  large-scale  fluctuations  in  w' .  At  higher  winds 
the  spatial  coherence  of  w'  is  destroyed  by  strong  turbulence.  Note  each  panel  in  hgure  4  is  sam¬ 
pled  at  the  sample  height  above  the  waveheld.  Also,  the  helds  are  made  dimensionless  by  fric¬ 
tion  velocity  u*  which  further  illustrates  the  strong  impact  of  the  wave  motions  on  the  winds  in 
the  surface  layer. 


In  figure  5  we  compare  vertical  profiles  of  the  mean  wind  speed  and  turbulence  variances  for  the 
different  simulations.  These  statistics  are  computed  by  averaging  in  computational  coordinates, 
i.e.,  across  horizontal  planes  at  constant  vertical  height  (.  Similar  to  our  previous  simulations 
we  find  that  the  wind  speed  and  turbulence  variances  depend  on  wave  age.  At  high  winds  as  the 
simulations  approach  wind- wave  equilibrium,  the  non-dimensional  wind  profile  (t/)/u*  smoothly 
approaches  the  variation  predicted  by  law-of-the-wall.  Significant  differences  are  observed  for 
the  cases  dominated  by  swell:  the  surface  layer  winds  are  accelerated  compared  to  rough  wall 
scaling.  As  suggested  by  the  flow  visualization,  the  turbulence  variances  respond  to  the  wave 
motion  in  dramatic  ways.  The  horizontal  and  vertical  variances  are  significantly  enhanced  by 
the  motion  of  the  wave  surface  in  the  low-wind  cases.  Even  though  the  turbulence  is  relatively 
weak  the  turbulence  variances  are  large  near  the  wave  surface  due  to  wave  pumping. 


IMPACT  /APPLICATIONS 


The  computational  tools  developed  and  the  database  of  numerical  solutions  generated  will  aide 
in  the  interpretation  of  the  observations  gathered  during  the  past  HIRES  field  campaign.  In  ad¬ 
dition  idealized  process  studies  performed  with  the  simulations  have  the  potential  to  improve 
parameterizations  of  surface  drag  under  high  wind  conditions  in  large  scale  models.  Also  future 
computations  will  utilize  wave  fields  measured  from  a  conventional  marine  X-Band  radar,  i.e., 
WAMOS  see  http://www.oceanwaves.org. 


TRANSITIONS  &  RELATED  PROJECTS 

The  current  work  is  a  collaborative  effort  between  NCAR,  numerous  university  investigators  and 
international  research  laboratories.  Also  the  present  work  has  links  to  the  ONR  DRI  on  the  im¬ 
pact  of  typhoons  in  the  Western  Pacific  Ocean  (ITOP). 
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Figure  1:  A  snapshot  of  the  waveheld  height  h(x,y,t )  that  is  imposed  at  the  bottom  of  the 
LES  code,  h  is  built  from  a  sum  of  linear  plane  waves,  and  waves  propagate  left  to  right  ac¬ 
cording  to  the  dispersion  relationship.  The  horizontal  grid  spacing  matches  the  LES,  i.e., 

Ax  =  Ay  =  2.34  m.  The  color  bar  is  in  units  of  meters. 


Figure  2:  An  instantaneous  x-z  slice  of  the  3D  time  varying  computational  mesh  in  the  lowest 
portion  of  the  PBL.  The  horizontal  gridlines  become  level  surfaces  at  about  100  m  above  the 
water.  Only  a  fraction  of  the  grid  is  displayed. 


Figure  3:  Snapshot  of  static  pressure  fluctuations  p'/p  in  an  x-z  plane  near  the  water  surface. 
The  upper  panel  is  a  swell  dominated  regime  with  wave  age  ~  4.8  while  the  lower  panel  is  a  case 
near  wind-wave  equilibrium  with  wave  age  ~  1.4.  The  wave  spectrum  is  a  Pierson-Moskowitz 
spectrum.  Notice  the  coherence  between  the  wave  held  and  the  pressure  huctuations  in  the  case 
with  swell.  The  color  bar  is  in  units  of  m  s~2  and  the  range  is  different  between  the  two  cases. 
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Figure  4:  Snapshot  of  resolved  vertical  velocity  fluctuations  w' /u*  in  a  wave  following  x-y 
plane  near  the  water  surface  (=2.5  m.  The  left  panel  is  a  swell  dominated  regime  with  wave 
age  ~  4.8  while  the  right  panel  is  a  case  near  wind-wave  equilibrium  with  wave  age  ~  1.4.  The 
wave  spectrum  at  the  bottom  of  the  PBLs  is  the  same.  Notice  the  range  of  the  color  bar  is  dif¬ 
ferent  between  the  two  cases.  The  (normalized)  fluctuations  in  the  wind-wave  equilibrium  case 
are  smaller. 


Figure  5:  Vertical  profiles  of  wind  speed  (left  panel)  and  turbulence  variances  (right  panels)  for 
different  values  of  wave  age  Cp/U\q.  Friction  velocity  u *  is  used  for  normalization.  The  dashed 
black  line  is  the  rough  wall  formula  U/u*  =  In (z/z0) /k,  where  n  =  0.4.  Temporal  and  spa¬ 

tial  averaging  is  used  to  make  the  statistics.  The  spatial  averaging  is  over  horizontal  planes  in 
computational  space,  i.e.,  at  constant  (. 


